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Abstract 

We emphasize intertwining relations as a universal tool in constructing one- 
dimensional quasi-exactly solvable operators and offer their possible generalization to 
the multidimensional case. Considered examples include all quasi-exactly solvable op- 
erators with invariant subspaces of monomials. We show that the simplest case of 
generalized intertwining relations allows to naturally relate quasi-exactly solvable op- 
erators with two invariant monomial subspaces to a nonlinear parasupersymmetry of 
second order. Quantum-mechanical systems with linear and nonlinear supersymmetry 
are discussed from the viewpoint of quasi-exact solvability. We construct such a general 
system with a cubic supersymmetry and argue that quantum-mechanical systems with 
nonlinear supersymmetry of fourth order and higher are generally not quasi-exactly 
solvable. Besides, we construct two examples of quasi-exactly solvable operators with 
invariant subspaces which cannot be reduced to monomial spaces. 
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1 Introduction 



In the eighties a new class of quantum-mechanical problems was discovered. Spectra of such 
systems can be partially found by an algebraic procedure, therefore they were called quasi- 
exactly solvable systems. In this sense these problems occupies an intermediate position 
between exactly solvable and unsolvable systems. There were elaborated several approaches 
to constructing quasi-exactly solvable quantal systems. Among them the Lie-algebraic con- 
struction P and analytic approach [2] are most common. In the Lie-algebraic approach a 
quasi-exactly solvable Hamiltonian is a polynomial of differential operators which form a 
finite-dimensional representation of a Lie algebra. This automatically guaranties the quasi- 
exact solvability of the Hamiltonian. In one-dimensional quantal problems the corresponding 
algebra is sl(2, R). Later on such ID systems were generalized to quasi-exactly solvable op- 
erators with invariant monomial subspaces |3] (for recent discussions see Refs. [3JE])- 

Recently it was revealed nontrivial relations between quasi-exactly solvable systems and 
nonlinear supersymmetry [HI 13 El El which is a generalization of the linear supersymmetric 
quantum mechanics [Tn | llT | IT2] to the case of nonlinear (polynomial) superalgebras [T8| I14j. 
The correspondence between the Lie-algebraic scheme based on the algebra sl(2, R) and su- 
persymmetric systems with nonlinear superalgebras was shown in Ref. [H] while supersym- 
metric systems corresponding to the quasi-exactly solvable operators with invariant mono- 
mial subspaces were constructed in Ref. |15j . 

Intertwining relations underline the standard realization of the quantum-mechanical su- 
persymmetry both linear and nonlinear. The only difference of the nonlinear case is that 
the corresponding intertwining differential operators are of higher orders. Intertwining rela- 
tions are corner stone of various approaches such as shape invariance to constructing exactly 
solvable systems |T2*1 ITS] . In this paper we emphasize the role of intertwining relations in 
constructing quasi-exactly solvable operators and offer a generalization of the construction 
to the multidimensional case. 

The paper is organized as follows. In the next section we discuss the relation of inter- 
twining relations to quasi-exact solvability of one-dimensional operators. Besides we offer a 
natural algebraic generalization of intertwining relations to the multidimensional case. In 
the section El we consider definition and some aspects of the one- dimensional quasi-exactly 
solvable operators. In the section 0] the results of the Lie-algebraic scheme based on the 
algebra sl(2, R) are reproduced in terms of the method based on the intertwining relations. 
In the section El we discuss application of the method to quasi-exactly solvable operators 
with invariant monomial spaces. An intimate relation of quasi-exactly solvable operators 
with two invariant monomial subspaces to a nonlinear parasupersymmetry of second order 
is outlined in the section El The section [7| is devoted to investigation of one- dimensional 
quantum-mechanical systems with linear and nonlinear supersymmetry from the viewpoint 
of quasi-exact solvability. In this section we construct a general quasi-exactly solvable oper- 
ator with three-dimensional invariant subspace, which is related to a cubic supersymmetry. 
Besides, using the method based on intertwining relations two new examples of quasi-exactly 
solvable operators are constructed. The corresponding invariant subspaces do not belong to 
spaces with a monomial basis but can be considered as their one-parametric deformation. 
The conclusion discusses and summarizes the obtained results. 
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2 Intertwining relations and quasi-exact solvability 



In a general sense a differential operator is called quasi-exactly solvable if it has a finite- 
dimensional invariant subspace. If this subspace is known explicitly the corresponding part 
of spectral problem can be solved algebraically. 

Let us consider the following intertwining relations 

AH = Hi A, A^Hi = H Q A\ (1) 

where H , Hi and A are some operators on a Hilbert space TL and f denotes an involution, 
e.g. a hermitian conjugation. As is well known the intertwining operator A relates the 
spectra of the operators Hi and H : 

IpE = Alp E 

for ipE ^ KerA = {ip e TC | Aip = 0}, where E stands for an eigenvalue. Besides, from the 
intertwining relations $1} it follows that if Ker A is not empty then it is a nontrivial invariant 
subspace of the operator H : 

Ker A Ker A 

This property is very interesting from the viewpoint of quasi-exact solvability. Indeed, if 
Ker/L is a finite-dimensional space then the operator H is quasi-exactly solvable. Hence, for 
a given operator A with a finite number of zero-modes the intertwining relations (0) provide 
a constructive way to build corresponding quasi-exactly solvable operator H . 

When the operators H Q and Hi are hermitian a superalgebra 1 can be associated with the 
intertwining relations (JTJ). Indeed, introducing the pair of odd variables 9, 9\ {9, 9^} = 1, 
one can construct the supercharges and super-Hamiltonian: 

Q = A9, Q f = A^9\ H = H Q 99 ] + H x 9 ] 9. (2) 

The supercharges commute with the super-Hamiltonian due to the intertwining relations (0) . 
Operators H, Q, generate a linear superalgebra fOJ when A is a linear differential 
operator. If A is a higher order differential operator the corresponding superalgebra is of 
nonlinear (polynomial) form [UJ [7| |H1 El HE] • 

The intertwining relations (0) are very useful for constructing one- dimensional quasi- 
exactly solvable operators since A being a one-dimensional differential operator always has a 
finite number of zero-modes. In multidimensional problems it is not the case and the scheme 
based on the intertwining relations should be modified in appropriate way. To this end we 
offer the following generalization of the intertwining relations (Q): 

M 

AiH = E H A v f o r i = 1, . . . , M. (3) 

i=i 

lr The operators also can be pseudo-hermitian and in this case one can associate a pseudo-supersymmetry 
Q21 E| with them. 



<r(H ) -^<t(HA i.e. 
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In this case if the operator- valued matrix Hij is "irreducible", i.e. it is not of a block- diagonal 
form, and Ha 7^ for all i, the space 

M 

F(H ) = f)KerA t ( 4 ) 

i=i 

is invariant subspace of the operator Hq. In general, depending on the particular form of the 
operator- valued matrix there can be some other invariant subspaces of H Q . For example, 
consider following two cases: 

A 1 H = H 1 A 1 , A 2 H = H 2 A 2 , (5) 

and 

A X H = H n Ai, A 2 H = H 21 Ai + H 22 A 2 . (6) 

In the former case if the intersection of Ker A{ is not empty there are three invariant subspaces 
of the operator H : Ker A4 \ (Ker A\ fl KerA 2 ), i — 1, 2 and Ker A\ D KerA 2 . In the latter 
there are two invariant subspaces: KerAi \ (Ker A\ fl Ker A2) and KerAi fl Ker^. More 
detailed analysis of the generalized intertwining relations Q in the context of invariant 
subspaces of the operator H is out of the scope of the present paper and will be given 
elsewhere. We just would like to add that although in the multidimensional case any space 
of zero-modes Ker Ai generally has infinite dimension the invariant subspace (j3J) can be 
finite-dimensional. Therefore the relations can be used for constructing multidimensional 
quasi-exactly solvable systems. 



3 Quasi-exactly solvable operators in one dimension 

Let us consider the most general differential operator of second order: 

H= -P( z )d 2 z +Q(z)d z + R(z), (7) 

where we assume that zEl = {z£M. \ P(z) > 0}, d z = d/dz, P(z), Q(z) and R(z) are 
some real functions. Formally, the operator H is hermitian with respect to the following 
inner product: 

where <fi(z) and ip{z) are real functions. The generalization to the supersymmetric case is 
straightforward. The wave functions become two component functions and the weight func- 
tion p(z) should be changed for a diagonal matrix defined from hermicity of the Hamiltonians 
Ho, H\. 
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We require for the Hamiltonian (J7J) to be a covariant operator with respect to coordinate 
transformations. Therefore, under the change of variable z — *> x(z) the functions P{z), Q{z) 
and R{z) transform as 

P(x(z)) = (d z x(z)) 2 P(z), 

Q{x{z)) = d z x{z)Q{z) - d 2 x(z)P(z), (9) 
R(x(z)) = R(z). 

When matrix elements of the Hamiltonian ((7|) are invariant under the coordinate transfor- 
mations: 

where we assume that the wave functions are normalized and transform as scalar fields, 
ip(x(z)) = ip(z), while the weight function p transforms according to its definition (JHJ). 

Besides, the system (J7J), (jHJ) has another local symmetry. The matrix elements are invari- 
ant under the so-called "gauge" transformations: 

$(z) = e a{z) i/j(z), H = e a{z) He~ a(z) . (10) 

In terms of the coefficient functions of the Hamiltonian it reads as 

P{z)=P{z), 

Q(z)=Q(z)+2d z a(z)P(z), (U) 
R(z) = R(z) - d z a(z)Q(z) + (d 2 z a(z) - (d z a(z)) 2 ) P(z). 

Form the definition (JHJ) one can deduce the transformation for p(z): 

p{z) = Xe- 2 ^ P (z). 

The constant factor A depends on the integration constant c in the definition of p(z) (jHJ) but 
it is irrelevant until we use normalized wave functions. 

The Hamiltonian ((7j) and the inner product (jHJ) can be represented in an explicitly co- 
variant form with respect to the transformations Q and To this end, consider the 
following operator: 

H = - -^V Z (^"V Z ) + V(z), (12) 

where the metric is given by the quadratic form ds 2 = g zz (z)dz ® dz, g = det||^|| = 
9zz , g zz gzz = 1 and V(z) is a scalar function. V z = d z — A z , where A z is component of 
the one-form A = A z (z)dz. The definitions of the metric and the one-form fix how their 
components change under the coordinate transformations while we assume that under the 
"gauge" transformations (fTT)|) the ID vector potential transforms as an Abelian connection 

A Z = A Z + d z a. (13) 

Therefore, the transformations (fTTj) are true gauge transformations and V 2 is the correspond- 
ing covariant derivative. Hence, the Hamiltonian (fT2*j) has an explicitly covariant form. 
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The appropriate invariant inner product has the form 

(0, V) = / <f>(z)iP(z)e- 2 f ZA ^gdz, (14) 

where the exponential factor provides the invariance of the inner product with respect to the 
gauge transformations (fTOjl. (fT5j) . 

The link between the covariant formulation (jEJ, |l4jl is given by the relations 

P(z) = g zz , Q(z) = 2g zz (A z + T z zz ) , R(z) = V(z) - A Z A Z + g zz V Z A Z1 

where Y z zz = d z log ^fg. So, there is a one-to-one correspondence between the coefficient 
functions of the Hamiltonian ((7j) and g zz , A z , V(z). 

Now let us discuss one-dimensional quasi-exactly solvable operators. If the operator ((Z|) 
is quasi-exactly solvable it has a finite-dimensional invariant subspace, say, 

T{H) = spanfACz), f 2 (z), f n (z)}, (15) 

where linearly independent real functions fi(z) form a basis of the subspace. One can think 
that the functions are ordered by number of their nodes. Hereinafter we will denote a finite- 
dimensional invariant subspace of an operator O as T(0). 

In principle one can extend this definition and treat an operator as quasi-exactly solv- 
able one if a part of its spectrum (such part can be infinite-dimensional) can be found 
algebraically. Here we discuss only quasi-exactly solvable operators with finite-dimensional 
invariant subspace. 

Generally, by a change of the coordinate and the gauge transformation (fTUJ) one can 
reduce the invariant subspace (fTKJ) to the following form 

JF(#) = span{l, z, f 3 (z), U(z), . . . , /„(*)}, (16) 

From this form it is obvious that the invariant subspaces of dimensionality 1 and 2 are special 
cases and require a special treatment. 

At the end of this section we would like to briefly discuss the general structure of quasi- 
exactly solvable operator. As was proved in Ref. ^H] every linear or nonlinear quasi-exactly 
solvable operator with a quite general finite-dimensional invariant subspace (fl3j) can be 
represented in the following form: 

n I 
i=l a=l 

where the e C w (M. n ), O a are arbitrary operators, A a are annihilating operators (annihila- 
tors) of the space (|T5|). i.e. A a [ip] = for all ip G JF(T 9e;s ), and Li form a basis of the space 
of affine annihilators, i.e. Li[ip] = const e M. for all ip G J~(T qe3 ). 

4 Intertwining relations and the $1(2, M)-based scheme 

In this section we consider application of the intertwining relations (0) to the well-known 
case of quasi-exactly solvable operators corresponding to the Lie-algebraic scheme based on 
the algebras 1(2, R). 



6 



The Hamiltonians H and Hi (0) are taken in the most general form 

Hi = - P^dl + Q<(z)0 a + R^z), (18) 

where Pi(z), Qi{z) and Ri(z) are arbitrary real functions. Here we treat the intertwining 
operator A as an annihilation operator of the corresponding invariant monomial subspace. 
Hence, the "bosonic" Hamiltonian Hq is a candidate for a quasi-exactly solvable operator 
while the operator H\ is its supersymmetric partner in the sense of the supersymmetric 
construction (JTJ), (J2J). 

In the case of the sl(2, R)-based scheme the invariant subspace can always be reduced to 
the form 

F n = span{l,z,z 2 ,...,z n - 1 }, (19) 

where n G N. It is the representation space of the algebra sl(2, M) with generators given by 
the differential operators 

71 — 1 

J+ = *d M -{n-l)z, J° n = zd z -—, J- = d z . (20) 

Here we assumed that n > 2 since the cases n — 1, 2 should be analyzed separately. 

In the standard approach any quasi-exactly solvable operator is taken as a polynomial 
of the generators (J22j) . However we will obtain the general form of such an operator using 
the strategy based on the intertwining relations This method is more general then the 
Lie-algebraic approach since, as we will see later, it is equally good for constructing any 
one-dimensional quasi-exactly solvable operators. 

It is very easy to see that an annihilator of the monomial subspace (|19|) can be represented 
in the following simple form: 

A n = d n z . (21) 

If we find the operators Ho, Hi and A n which obey the intertwining relations then Hq is a 
quasi-exactly solvable operator with the finite-dimensional invariant subspace KerA n = T n . 
The question on the normalizability of such states we leave aside in this paper. 

Substituting the operators (fTSjl and (J2T|) into the intertwining relations (JTJ) one gets a 
system of differential equations for the functions Pi{z), Qi{z), and Ri(z), i = 1, 2. This 
system can be divided into two parts. The fist one can be solved algebraically and allows to 
express the coefficient functions of Hi in terms of those of Hq: 

P 1 (z)=P (z), 
Q 1 (z)=Q {z)-nP '{z), 

71 ( T) — 1 I 

Ri(z) = R (z) + nQ' (z) - A__d P » {z) . 

The second part is an overdetermined system of differential equations for the functions Pq(z), 
Qo(z) and R (z): 

m(m - l)P (n_m+2) (^) -m(n-m + 2)Q { n - m+1 \z) - (n - m + l)(n - m + 2)R { n ' m \z) = 0, 
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where m = 0, . . . , n — 1. The equations with m G {n — 3, n — 2, n — 1} are independent and 
give the following 9-parametric solution: 

4 

p ( Z ) = Y,p^ m ^ 

Qo(z) = ^q m z m + T ^P[ ) {z) 1 (22) 

(n—l)(n — 2) ... . n — 1 , . . 
Ro{z) = K - >j L PS(z)-— Q' (z)+v, 

where p m , q m , and v G R. This solution satisfies the complete overdetermined system of 
differential equations. The second order operator Hq with the coefficient functions (}2*2*j) has 
the same form as in the Lie-algebraic scheme based on the algebra sl(2, R). 

It is worth noting that the procedure of solving such a system of differential equations is 
so simple only in the case (|19j) . 



5 Generalized quasi-exactly solvable systems with in- 
variant monomial subspaces 

A generalization of one-dimensional quasi-exactly solvable systems corresponding to the 
st(2, R)-based approach was investigated for the cases of invariant subspaces with a mono- 
mial basis [3] (for more recent discussion, see Refs. jUE])- In this section we demonstrate 
that the approach to quasi-exact solvability based on the intertwining relations (JTJ) can be 
applied to such quasi-exactly solvable systems as well. 

Now let us discuss how one can deal with a general monomial invariant subspace using 
the intertwining relations (0). An n-dimensional space with monomial basis functions can 
be represented as 

F = span{z mi , z m2 ,..., z mn }. (23) 

One can consider the parameters m, as natural numbers but for the method based on the 
intertwining relations it does not matter and rrii can be treated as real numbers. 
The annihilator of the space ()23|) can be taken in the form 

n 

A n = Y[{zd z - mi ). (24) 

i=l 

The intertwining relations (0) with the intertwining operator A n lead to the overdetermined 
system of differential equations on coefficient functions of the operators Hq and Hi. The 
particular solution of the system depends on the parameters m ; and order of the operators 
Hq and Hi. A solution always exists since H taken as a polynomial of the operator zd z 
leaves the subspace (|23|) invariant and the form of coefficient functions of the operator Hi 
can be found from the intertwining relations algebraically. We are interested in nontrivial 
solutions for the differential operators Hq and Hi of second order. 
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In terms of the new variable £ = log z one has zd z = and the annihilator (|24j) becomes a 
differential operator with constant coefficients. As a consequence, the intertwining relations 
(J!Q) lead to a system of differential equations with constant coefficients for the coefficient 
functions of the Hamiltonians Hq and H\. A system of this kind can be solved by a purely 
algebraic procedure. Therefore the resulting Hamiltonian Hq is a quasi-exactly solvable 
operator with the invariant subspace ([23)1 . 

Here we will not elaborate this procedure in any detail. Linear differential operators with 
finite-dimensional invariant subspaces of monomials were investigated in Ref. [3] 2 , where a 
classification of second order operators with invariant subspaces of the form (|23|) was given. 

We give here the classification of such finite-dimensional invariant subspaces with basis of 
monomial functions without any explicit form of the corresponding operators. The form of 
the operators can be found in Refs. ISHHEj- Invariant subspaces of monomials for nontrivial 
second order differential operators 3 can be divided into two classes. The fist one comprises 
the following spaces with a monomial basis of dimension 3 and 4: 

T z = spanjl, z, z m }, JF 4 = spanjl, z, z m , z 2m ). 

The parameter m can be treated as real number. The second class is represented by series 
of spaces of monomials. The first series is that corresponding to the Lie-algebraic approach 
(|19p . The next series is 

T = spanjl, z, z 2 , z n ~ 2 , z n }, (25) 
where n G Z. The last series is a direct sum of two subspaces: 

T = spanjl, z, z 2 , . . . , z n ~ 1 } © span{z m , z m+ \ z m+2 , z™^- 1 } = T n © z m F k , (26) 

where n, k G N and m G R. 

Up to a change of variable and gauge transformations these spaces exhaust the list of 
invariant monomial subspaces of quasi-exactly solvable differential operators of second order. 

All these cases can be reproduced in terms of the usual intertwining relations ((TJ. This 
method allows to automatically construct the associated supersymmetric systems with non- 
linear superalgebras (details on the construction can be found in Ref. jlH])- in the next 
section we will show that using the simplest generalized intertwining relations (0) a non- 
linear parasupersymmetry can be naturally associated with quasi-exactly solvable operators 
invariant on the monomial spaces (|26|). 

6 Systems with two invariant subspaces and nonlinear 
parasupersymmetry 

In quasi-exactly solvable systems with two invariant subspaces f!26|) one can associate an 
annihilating operator with each subspace. Consider the following two operators: 

A 1 = d n z} A 2 = d k z z~ m . (27) 

2 More recent discussion on the subject can be found in Refs. 0J|S]. 
3 Here the trivial operator is a polynomial in zd z . 
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Their spaces of zero-modes are 



Ker/l! = T n , Ker A 2 = z m T k . (28) 

If the Hamiltonian H is a second order differential operator and obeys the generalized 
intertwining relations (jSJ for some operators H% and H 2 then it is a quasi-exactly solvable 
operator with the invariant subspace (j2E| 

T{H) = Ker Ax © Ker A 2 =J r n ® z m T k . 

From the first relation in (j5J) it follows that the Hamiltonian H is a quasi-exactly solvable 
operator corresponding to the sl(2, M)-based scheme. The second relation from (J5J) can be 
considered as a constraint on this quasi-exactly solvable operator. Thus, without explicit 
calculations we came to the same conclusion as made in Ref. We can even say something 
more about this system. It is known that the generalized intertwining relations (0) underlay 
a parasupersymmetric system [20] • Indeed, from the relations (0) one concludes 



AlA^P^Ho), A i A\ = P i (H, 



1 / ■ 



where i — 1, 2, P\(.) is a polynomial of order n and P 2 (-) is a polynomial of order k. Here the 
hermitian conjugation is understood with respect to the quite straightforward generalization 
of the inner product (JHJ) to the case of 3-component wave functions and p changed for a 
3x3 diagonal matrix defined from hermicity of the Hamiltonians H , Hi, and H 2 . If one 
introduces the new operators 

ai = Ax, a 2 = A\ (29) 

the spectra of the Hamiltonians Ho, Hi and H 2 are related as follows 

<t{H x ) A a(H ) A a(H 2 ), <r(H x ) ^ a(H ) a(H 2 ). 
Let us introduce the following matrix operators 

'H 2 

H = I H I , Q=\a 2 1. (30) 




H h 

These operators form the nonlinear algebra given by 

[H, Q] = 0, Q 2 Q ] + QQ ] Q + Q ] Q 2 = Q(P 1 (H) + P 2 (H)), Q 3 = 0, (31) 

and the hermitian conjugated relations. It is easy to see that this algebra is a nonlinear 
version of the linear parasupersymmetry of order p = 2 [20] • This allows to conclude that the 
one-dimensional quasi-exactly solvable systems with two invariant subspaces are naturally 
related not only to the nonlinear supersymmetry but to the nonlinear parasupersymmetry 
as well. 
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7 Supersymmetric quantum mechanics from viewpoint 
of quasi-exact solvability 

As it was discussed in Section |21 the intertwining relations (JTJ) underline the general con- 
struction of one-dimensional supersymmetric quantum mechanics. If the operator A is a 
linear operator one arrives to the usual supersymmetric quantum mechanics [HH ITT] while 
the construction with A to be a higher order differential operator corresponds to a nonlinear 
generalization of the supersymmetry with a polynomial superalgebra [SJ El EE] • Owning 
to the discussed above correspondence between the intertwining relations (JTJ) and quasi- 
exact solvability of one-dimensional systems we can consider the supersymmetric quantum 
mechanics in this context. 

In a general quantum mechanical system with a linear supersymmetry the "bosonic" 
Hamiltonian can be treated as a quasi-exactly solvable operator with a one-dimensional 
invariant subspace. Since the sequence of eigenfunctions of the Hamiltonian can always be 
reduced to the form (jl6j) without loss of generality we can take the invariant subspace in the 
form T{H) = {1}. It is easy to see that if the second order differential operator is invariant 
on this space then it has the following form: 

H= - P(z)d 2 z +Q(z)d z + const. (32) 

After changing the variable 



x 



dz 



and introducing the notation 



W{x) 



P'{z) + 2Q{z) 



z=z(x) 



A^Pjz) 

the Hamiltonian (]3*2"|) acquires the standard form 

H = -dl + W(x) 2 - W\x) + const. 



One can note that while the operator (|32j) depends on two arbitrary functions the resulting 
operator of the Schrodinger form counts on just one function, the superpotential. 

It is interesting to note that in this case the set of annihilating operators is infinite- 
dimensional. Its basis can be chosen in the form L n = z n+l d z , n G Z, hence, these operators 
form the Virasoro algebra without central extension. The operator (|32|) can be treated as 
a combination of quadratic and linear terms in the generators of the Virasoro algebra. One 
can conclude that this fact is responsible for the presence of the functional arbitrariness in 
the supersymmetric Hamiltonian. 

Now consider the general case of a second order quasi-exactly solvable operator with two- 
dimensional invariant subspace. Without lost of generality one can deal with the following 
subspace 

T{H) = span {1, z} . (33) 
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The corresponding quasi-exactly solvable Hamiltonian can be written as 

3 

H= - P(z)d 2 z + C aJi/2 + const, (34) 
0=1 

where C a are constant and J®, 2 are linear differential operators fl22(l forming a representation 
of the algebra sl(2, M) on the two-dimensional space (J33j) . 

The structure of the operator (|3"4*|) with the two-dimensional invariant subspace is dif- 
ferent as against the previous case. The second term in ()34j) is a linear combination of 
affine annihilators of the invariant subspace while the first one can be treated as a linear 
combination of annihilators for that space. 

The annihilators of the space ()33p are any differential operators of order 2 and higher 
without linear in derivative and multiplicative terms. Obviously, they form an infinite- 
dimensional set. The basis of this set can be taken in the form 

L n = z n+2 d 2 z , neZ. (35) 

These operators together with the affine annihilators Jy 2 generate the following nonlinear 
infinite-dimensional algebra: 

\_L m , L n ] = (m — n)L m+n o J{y 2 > [^1/2 > ^m] = m L m , [jjy 2 ,L m ] = (m =p 2)L m ±i, (36) 

where o stands for anticommutator. Therefore, one can conclude that the annihilators (|35j) 
together with the affine annihilators form the nonlinear (quadratic) generalization of the 
Virasoro algebra. 

Using the usual intertwining relations (£[]) with the annihilator A = d 2 z one can reproduce 
the form of the Hamiltonian (}3"4*|) and construct its superpartner. That would complete the 
construction of the supersymmetric system with a quadratic superalgebra [T5] . 

Now let us investigate the case of a three-dimensional invariant subspace, which can 
always be reduced to the following form: 

^(if) = span {1, *,/(*)}. (37) 

Because of the presence of the arbitrary real function f(z) this case looks different in com- 
parison with the two previous cases. As we shall see it is also different from the viewpoint 
of quasi-exact solvability of second order operators. 

Using the usual intertwining relations (|1J) with the annihilator of the subspace (J37|) 

A = dl-{\ogf{x))'d 2 z (38) 

one can derive the form of corresponding Hamiltonian: 

3 

H=Y f f i {z)a ij K j , (39) 
i,i=l 

where the functions 

fi(z)=l, Mz) = z, h{*) = f{z) (40) 
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form the basis of the invariant subspace ()37|) and are real numbers. The operators Ki 
form a basis of affine annihilators of the subspace (|37j) : 

Kl = l-zK 2 - f(z)K 3 , K 2 = d z - j^fl K z = f\z)- l dl 

They act on the basis ()40|) as follows: 

Ki[fj] = 5ij. 

From the viewpoint of the general structure (|17|) of quasi-exactly solvable operators the 
Hamiltonian ()39|) is just a linear combination of affine annihilators. Any annihilator of the 
invariant subspace ()37|) is a differential operator of order not less then 3. Therefore, the 
Hamiltonian cannot contain any of them since we are restricted by the class of second order 
differential operators. 

The form of the supersymmetric partner of the Hamiltonian (jHUjl can be derived from the 
corresponding intertwining relations. We do not need to know it to obtain the form of the 
superalgebra they generate together with the supercharges (see (J2J)). The conservation of 
supercharges follows from the intertwining relations. The anticommutator of supercharges 
is proportional to a cubic polynomial in the super-Hamiltonian: 

{Q, Qt} = Det (El - A)\ E ^ k = H 3 + . . . , 

where elements of the matrix A — II defined in (J39j) . This follows form the fact that 

the finite-dimensional invariant subspace (}3*Tj) of the "bosonic" Hamiltonian (|39J) corresponds 
to the zero-mode space of the supercharge Q. 

Construction of the cubic supersymmetry with Hamiltonians in Schrodinger form was 
presented in Ref. |2T]. Other approaches to general quasi-exactly solvable operators with 
three states can be found in Ref. 

Now we would like to discuss the general case of one-dimensional quasi-exactly solvable 
differential operators of second order. If the operator (JJJ) has a n-dimensional invariant sub- 
space (fT5j) its coefficient functions and the basis functions f%{z) obey the system of differential 
equations 

n 

- P(z)f7(z) + Q(z)fl(z) + R(z)f t (z) = J2 C afj(z), (41) 

z = 1, ... ,n, for some real constants CV,-. If n < 3 the equations can be treated as a 
system of algebraic equations for the variables P(z), Q{z), and R(z). In this manner the 
Hamiltonians and (jH9*|) can be reproduced. Beginning with n = 4 the situation 

changes drastically. Indeed, in the case n = 4 the system allows not only to fix the coefficient 
functions of the Hamiltonian but it also leads to a nonlinear differential equation for the 
functions fi(z). 4 

Besides, from the system (J41j) one can easily infer that in the cases n — 1, 2 the func- 
tional arbitrariness of the Hamiltonians is related to the degeneracy of the corresponding 
algebraic system. The system (}4"Tj) also allows to conclude that the nature of the functional 

4 Or for the functions fo(z) and f±(z) in the case of the reduced invariant subspace (|16fl . 
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arbitrariness in the Hamiltonian (J3U|) is different as against the cases ()32|) and (|34p. For 
n = 1 and n = 2 the functional arbitrariness is related to the degeneracy of the system of 
algebraic equations on the coefficient functions while for n = 3 it is related to the presence 
of the arbitrary function in the basis (jl6j) . 

The Hamiltonian restricted to the invariant subspace is represented by the matrix dj. 
Although it is real valued but nonsymmetric matrix, hence, its eigenvalues can be complex 
numbers. Therefore, if CV, is nonsymmetric matrix one has to assume that the functions 
forming the basic of the subspace do not have finite norm with respect to the inner product 
(JB1) since the Hamiltonian has to be a hermitian operator. General structure of the matrix 
Cij was studied in Ref. [THj in the context of polynomial superalgebras. 

In the general context the most nontrivial problem is to find and classify all functional 
finite-dimensional spaces which can serve as invariant subspaces of differential operators 
of the general form (J7|) and cannot be reduced to the monomial form using changes of 
the variable and gauge transformations. Such spaces do exist. For example, consider the 
following one-parametric deformation of the series (|25|): 

T = spanjl, z, z\ . . . , z n ~ 3 , z n ~ 2 + -^—z n ~\ z n \ (42) 

n — 1 

with n G N and a G 1R, which cannot be reduced to the monomial form for a / 0. The 
corresponding annihilator of the subspace can be represented in the form 

A n = ((az + 2)d z - a) (zd z - 2)d n z ~ 2 . 

Application of the method of intertwining relations leads to a 6-parametric quasi-exactly 
solvable Hamiltonian which is particular case of the 5l(2, M)-based Hamiltonian with n — > 
n — 2 and the following constraints on the parameters: 

n + 2 2 n + 4 

Va = 0, g = —z— Pi ~ apo, qi= (n - l)p 3 + np 2 + ap 1 - a 2 p , q 2 = ~^—P3- 

z a 2 

As another example of invariant subspaces of monomials consider the following space: 

T = spanjl, z, z 3 , z 4 + Qaz 2 }. 

Using the method of intertwining relations one calculates the 7-parametric quasi-exactly 
solvable Hamiltonian given by the coefficient functions 

P(z) = p 5 z 5 + p 3 z 3 + p 2 z 2 + piz + p , 

Q{z) = 6p 5 z 4 + 4 (p 3 + 4ap 5 ) z 2 + (5p 2 - a _1 po) z ~ 8a 2 p 5 + 4ap 3 + 2p u 
R(z) = - 12p 5 z 3 - 4 (13ap 5 + Pa) z + r . 

More detailed discussion of quasi-exactly solvable operators with non-monomial invariant 
subspaces will be presented elsewhere. 
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8 Conclusion 



In this paper we emphasized the method based on the intertwining relations as a general ap- 
proach to constructing one- dimensional quasi-exactly solvable operators. We demonstrated 
the universality of this method, which is based on the algebraic nature of the intertwining 
relations. The discussed examples cover all quasi-exactly solvable operators with invariant 
monomial subspaces. It can be equally applied to quasi-exactly solvable operators with 
non-monomial invariant subspaces. Using the method we constructed two examples of such 
operators with invariant subspaces which cannot be reduced to a monomial space. 

The intertwining relations themselves can be applied to the one-dimensional case only. We 
offered natural algebraic generalization that can be applied to the multidimensional case as 
well. We showed that the simplest form of the generalization allows to relate the quasi-exactly 
solvable systems with two invariant monomial subspaces to a nonlinear parasupersymmetry 
of second order. The application of the generalized intertwining relations to constructing 
quasi-exactly solvable operators in the multidimensional case will be presented in future 
publications. 

Besides we investigated quantum-mechanical systems both with linear and nonlinear su- 
persymmetry from the standpoint of quasi-exact solvability. We demonstrated that in the 
cases of linear and quadratic supersymmetry the functional arbitrariness of the Hamilto- 
nian is related to an infinite-dimensional algebra of annihilators of corresponding invari- 
ant subspaces. It is the Virasoro algebra without central extension in the first case and a 
quadratic deformation of the Virasoro algebra in the second. Also we constructed a gen- 
eral quantum-mechanical system associated with a cubic supersymmetry and argued that 
quantum-mechanical systems with nonlinear supersymmetry of order more then three are 
not quasi-exactly solvable in general. 
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